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Abstract 

We present a density-functional study of a binary phase-separating mixture of soft core particles 
immersed in a random matrix of quenched soft core particles of larger size. This is a model for 
a binary polymer mixture immersed in a crosslinked rigid polymer network. Using the replica 
'trick' for quenched-annealed mixtures we derive an explicit density functional theory that treats 
the quenched species on the level of its one-body density distribution. The relation to a set of 
effective external potentials acting on the annealed components is discussed. We relate matrix- 
induced condensation in bulk to the behaviour of the mixture around a single large particle. The 
interfacial properties of the binary mixture at a surface of the quenched matrix display a rich 
interplay between capillary condensation inside the bulk matrix and wetting phenomena at the 
matrix surface. 
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I. INTRODUCTION 



Spatial confinement can have a dramatic effect on the properties of a fluid. The com- 
bined effects of finite pore size and substrate fluid interactions can alter the location of phase 
boundaries from what is found for the bulk adsorbate. Many studies have focused on the 
equilibrium structure and phase behaviour of fluids confined in model pores or capillaries and 
phenomena such as capillary condensation (the shifted gas-liquid transition) and capillary 
freezing are rather well-understood for a one-component fluid confined in a single, infinitely 
long slit or cylinder . The kinetics of nucleation accompanying capillary condensation 
has also been investigated . When the adsorbate consists of more than one species, capil- 
ary induced fluid-fluid demixing can occur in addition to capillary condensation - see e.g. 
^| for recent work. However, most porous media exhibit complicated, interconnected pore 
shapes P] and cannot be modelled as idealised capillaries. Rather, a random porous network 
structure provides a more realistic starting model. In recent years the direct treatment of 
fluids adsorbed random pore structures has been achieved within the framework of lattice 
fluids yielding significant new insight into adsorption phenomena 

An alternative way of modelling porous structures is to consider pore structures generated 
by quenched configurations of model fluids. These configurations carry a statistical weight 
that is determined solely by the Hamiltonian of the (pure) model fluid. The adsorbate fluid 
is then brought into contact with the quenched particles (referred to as the matrix), and is 
"annealed" in the sense that it is under the influence of the matrix as an external potential. 
In addition to the ensemble average, an average over all realisations of the matrix must be 
carried out. Frequently such models are tackled using liqu id state integral equation theory, 
relying on the replica Ornstein-Zernike equations [llj, [l2| - an approach that works on the 
level of two-body correlation functions. 



I '., [ljj], operating on 1 lie one 



Recently a classical density functional theory (DFT) 
body level, was shown to be able to treat consistently quenched-annealed mixtures. This 
quenched-annealed DFT (or replica DFT) treats both the annealed and the quenched 

species on the level of their one-body density distributions. Hence the complicated external 
potential that the quenched particles exert on the fluid never enters explicitly into the theo- 
retical framework. This approach constitutes an enormous simplification as far as practical 
computations are concerned. Investigations were carried out for hard core models and for 
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the quenched-annealed Asakura-Oosawa-Vrij model of colloid-polymer mixtures |l7j ]. 

The Gaussian core model (GCM) [18] belongs to a different class of system, characterised 
by fully penetrable pair interactions, i.e. at zero spatial separation between the centres of a 
pair of particles the potential energy is finite and is comparable in magnitude to the thermal 



energy. The fact that a very simp 



fluid densities 
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e density functional approximation is accurate at high 
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and much is known about the 
phenomena, in this model fluid 
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231 ] has facilitated a number of detailed investigations 
k and interfaci al p roperties as well as subtle solvation 
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In the present work, we are concerned with a quenched-annealed GCM modelling a binary 
polymer solution as an adsorbate in a quenched (crosslinked) rigid polymer network - i.e. a 
mobile binary mixture adsorbed in a model gel. Quenching components of a GCM mixture 
results in soft confinement of the annealed species, as opposed to the hard confinement 
considered with replica DFT in previous studies. The particularly simple (RPA-like) form 
of the excess free energy functional enables us to elucidate the structure of the replica DFT. 
We can relate the effect of the quenched component directly to the influence of an effective 
external potential acting on the fluid particles; this is shown to take the form of a convolution 
of the matrix one-body density profile and the matrix-fluid (pair) interaction potential. 

The paper is organised as follows. In Sec. |H] we give a brief overview of the statistical 
mechanics of quenched-annealed mixtures. In Sec. IHII the DFT for the quenched-annealed 
GCM is presented. The case of matrices with uniform density is covered in Sec. IIVI We 
present results for bulk properties of a particular mixture in SecEJand treat inhomogeneities 
induced by the surface of a matrix in Sec. IVII Discussion and conclusions are presented in 
Sec.lVTfl 



II. QUENCHED ANNEALED SYSTEMS 

We consider a system composed of three different species of particles. One of the species, 
labelled 0, is quenched, forming a matrix into which the remaining two species, labelled 1 
and 2, are annealed, i.e. the particles of species form an immobile network in which the 
other species are free to move. The correlations between the quenched matrix particles are 
independent of the adsorbate particles and are determined solely by interactions between 
matrix particles. We denote particle numbers by N a , where a = 0, 1, 2 labels the different 
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species, and denote the sets of position coordinates of particles of species a by = {r ija , i = 
1, . . . , N a }. The total external potential energy of species a, the total potential due to pair 
interactions between particles of species a, and the total potential due to pair interactions 
between particles of (unlike) species a and b, are defined, respectively, by 

N a 

K(r^) = ^K Cxt (r,), (1) 

N a N a 

<Mr" a ) = E E ^(l^-^l)' (2) 

i=l j=i+l 

N a N b 

ma^) = EE^-^)> a ^ b ( 3 ) 

i=l j=l 

where V a ext (rj) is the external potential acting on particle % of species a and v ab (\ri — rj\) is 
the pair potential between particle i of species a and particle j of species b. It is important 
to note that for quenched-annealed systems the potentials $oi and $02, exerted by the 
quenched particles on the annealed matrix particles, never enter into the description of any 
properties of the matrix (i.e. species 0) alone; the latter originate solely from V^r^ ) and 
$00 (r N °). 

Hence, the partition function for the quenched matrix particles is 

HoG^o) = Tr No exp[-/?($ 00 + V )\ (4) 
where /x is the chemical potential for the quenched particles and the classical trace is 

00 N (■ 

Tr ^(-)= Eli/*-). ( 5 ) 

where z a = A~ 3 exp(/3/x a ) is the fugacity and A a is the thermal de Broglie wavelength for 
species a, and (3 = l/fc^T, where ks is the Boltzmann constant. Strictly, S = S (/i ,T, V) 
is a function of the chemical potential /j,q, the temperature, T, and the volume of the system, 
V. However, we will suppress the dependence on T and V for notational convenience. The 
grand potential for the quenched particles is then given by 

^o(/V) = -k B T\nE (ii ). (6) 

For a given configuration of matrix particles, r^ , the grand potential for the annealed 
particles, adsorbed in the matrix, is: 

E 12 (r No , f i 1 , f i 2 ) = Tt Ni Tt N2 exp[-/3(Vi + V 2 + $ 12 + $ i + $02)], (7) 
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and the grand potential for this configuration is 



612(1^1,^2) = -k B T\nE 12 (r N ", f i 1 , f i 2 ) 



(8) 



We now average this grand potential over all configurations of quenched matrix particles: 



^12(^0,^1,^2) = H 1 Tr A r o ni 2 (r Aro ,/ii,/! 2 )exp[-/3($oo + V )]. 



(9) 



Q12 is a key quantity: it determines many of the properties of the fluid in the matrix. We 
now apply the replica 'trick' that proves to be a powerful idea to actually calculate f2 12 

The idea is to introduce s 'replicas', or copies, of the sets of particles of species 1 and 2 
into the matrix of quenched particles. Particles do not interact with particles of a different 
replica, but only with particles of the same replica as well as with the matrix particles. The 
partition sum for this (2s + l)-component mixture is: 



S(//o, HuH2\ s) = Tr No exp[-/3($ o + V )} {Tr^Tr^ exp[-/?(^ + V 2 + $ 12 + $ 01 + $ 02 )]} s 

(10) 

and the corresponding grand potential follows as 



(11) 



tt(no, Hi, fa] s) = -k B TlnE(iJ l0 ,Hi,H2;s). 



Performing an analytic continuation in s, and noting that lim s ^ dx s /ds = lnx, we find jl 
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^12(^0, Hu H2) = hm —Q(ho, Hu H2] s) 
s-»0 os 



(12) 



Thus we have constructed a means of tackling quenched-annealed systems by considering 
the corresponding replicated (and fully annealed) mixture. At this point we can make 
connection with DFT, which provides a powerful approach to calculating the grand potential 
for inhomogeneous systems. In DFT t he g rand potential is expressed as a functional of the 
fluid one-body density profiles, p a (r) llJ]. The grand potential in (fTTjl is obtained by 



considering a variational gr 'and potential functional, fi var ([{p a }], {Ha}] s), which is minimised 
with respect to the density profiles 



6n™(i{ Pa }],{tA a y,s) 

<*Pa(r) 



(13) 
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Inserting the density profiles that are the solutions to (JT3jl into Q vaT yields as the minimal 
value the true grand potential, fi(/i , pi, P2; s). Performing an analytic continuation in s, 
and Taylor expanding Q var around s = 0, gives 

fi var ([{pj], W; s) = fir(N^o) + *ng([{Pa}], W) + o(s 2 ). (u) 

^o ar ([Po] ; A*o) is simply the grand potential for the pure system of matrix particles; further- 
more ^([{ Pa }] 7 {/i a }) = lim^ ^ w ([{p a }],{/i a }; S )/fe. The key point is that Q\f in 
f!14|) corresponds to an average over the different matrix configurations and is hence only a 
functional of the matrix density profile Po(r), rather than depending on particular matrix 
configurations r^ . Explicitly, we can write the functionals on the right hand side of (|T4*j) as 

^([PoWUo) = n d [Po(r)] +-F CX [Po(r)] - J drp (r)[/i - V^{v)l (15) 

where 

r a d [p a (r)} = k B T f drp a (r)[ln(A^ a (r)) - 1] (16) 

is the Helmholtz free energy of the ideal gas and J r Q X [po(r)] is the excess contribution to the 
total Helmholtz free energy arising from interactions between the matrix particles. Similarly, 

^([p ,Pi,p 2 ],/ii,^ 

-J2 I drp a (r)K-K cx (r)], (17) 

a=l J 

where JF^ d [p a (r)], a = 1,2 is given by Eq. (fT7)|) and [po, pi, P2] is the contribution to 
the Helmholtz free energy from interactions between the particles. This quantity includes 
interactions between the annealed particles and the quenched matrix particles. In general, 
JF°2 is an unknown functional. However, depending on the form of the interactions between 
the particles, there exist a number of accurate approximate functionals '141. 
The variational principle (JEJ), when applied to ((HJ), becomes j 1 -"31 11 6j: 

^r(boUo) = 0j (18) 



<Vo(r) 
5^| r ([po,pi,p 2 ],Pi,P2) 



0, a = 1,2. (19) 



^Pa(r) 

One therefore uses (fTHj). together with (JT3J) and a suitable approximation for jFQ X [p (r)], 
to first calculate the average one body density for the quenched matrix particles, po(r). 



Then one uses this density profile as input to calculate the density profiles pi(r) and p2(r) 
of the annealed species using Eqs. (fTTjl and (fT9*|). again with a suitable approximation for 
jFg [{p a (r)}]. In the next section we shall describe specific approximations for particles that 
interact via repulsive Gaussian pair potentials. 



III. DFT FOR THE QUENCHED ANNEALED GCM 



The interactions between macromolecules such as polymers or dendrimers in solution 
can be modelled by means of an effective pair potential between the centres of mass of 
the macromolecules. This implies treating each macromolecule as a (soft) 'particle'. For a 
good solvent the effective p air p otential is purely repulsive and is well-approximated by a 
repulsive Gaussian form [3,130,0, 3 [3 Q- [3 • I n the GCM, the potential between 
particles of species a and b is given by: 



v a b{r) = e ab exp(-r 2 /R 



ab) 



(20) 



where e ab is the energy penalty (due to entropic effects on the segment level) for complete 
overlap of the centres of mass. Typically e ab ~ 2ksT and R aa ~ Rg,ai the radius of gyration 
of particles of the a th species. The pair potential parameters for the potentials between 
unlike-species, as obtained from simulations of model polymer solutions |29J], are consistent 
with the (non-additive) mixing rules, 

and 



-,( R U + R 22) 



(21) 



ei2 < en 



(-22- 



For the GCM the mean-field excess Helmholtz free energy functional, 
F CX [{Pa(r)}] = / dr / drV a (r)p b (r> a6 (|r-r'|), 



(22) 



(23) 
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23| when each particle inter- 



is known to be accurate at high densities 
acts with a large number of neighbouring particles - the classic mean field situation. The 
functional (}23|) generates the RPA closure for the pair direct correlation functions upon 
differentiation llj, Il8| since 



5p a (r)5p b (r') 



-(3v ab (\r-r'\). 



(24) 



Here we focus on the case of a ternary mixture where one species of GCM particles is 
quenched (treating cases with larger numbers of both annealed and quenched components 
requires straightforward generalisations). The quenched particles have a one-body density 
distribution, po(r), which is obtained by minimising (|15|) using the one component (0) ver- 
sion of (|23|) . The density distributions of the remaining (annealed) species of particles are 
obtained by minimising the grand potential functional: 

^2 r KPa(r)}] = ^ d [ Pl (r)]+^ d [p 2 (r)] 

+ \ ^2 / dr / dr'p n (r)p m (r')v nm (\r - r'\) 

n,m=l 

+ / dr / drVn(r)p (r>o„(|r-r'|) 

n=l ^ J 

-J2 [ drp„(r)[p n -\/rW]. (25) 

71=1 "* 

Expressing the density distribution of the quenched particles as po( r ) — Po + Ap (r), i.e. as 
a modulation Apo(r) around a constant density pg, ()25j) can be written as: 

U£[{pa{r)}] = ^ d [ Pl (r)] +^ d [p 2 (r)] +^ cx [{p n (r)}] 

-J2 I drp„(r)[pf -Vf(v)\ (26) 

n=l ^ 

where the effective chemical potentials are 

P,f = Pn-p^0n(0), (27) 

and t>on(0) = / drt> 0n (r) = 7T 3 ' 2 e n? i?^ 9 is the A; — > limit of the Fourier transform of the pair 
potential vo n (r). The effective external potentials are 

Vf(r) = K ext M + J dr'Ap (r> 0n (|r - r'|). (28) 

This analysis demonstrates that the quenched matrix has a very similar effect to that of an 
external potential acting on the mobile components. However, because of the convolution 
form in Eq. (|2*%|) . the matrix effective potentials V^ ff (r) are quite different from the external 
potentials that are typically tackled within DFT, such as the hard potentials modelling 
container walls. 
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IV. A MATRIX OF UNIFORM DENSITY 



As a prerequisite for our subsequent interface study we consider the case where the matrix 
density profile is constant: Po( r ) = Po (i- e - ^Po( r ) = 0). In the particular case where the 
external potentials also vanish, the one body densities of the annealed species are constants: 
p n (r) = , n = 1, 2. As convenient variables we use the total (bulk) density of the annealed 
species, p b = p\+p 2 , and the relative concentration, x = p\j p b . From Eq. (j2Mthe Helmholtz 
free energy per particle, f± 2 = F/N, for the annealed GCM fluid is simply [21 1: 

fn(p b ,x) = fH(p b ,x) + \p b V {x) + p b Q(x), (29) 

where is the ideal gas contribution, 

fH(p\x) = (1 -x)ln(l -x) +x\yix 

+ (1 - x) ln(A^ 6 ) + x \n(A 3 2 p b ) - 1, (30) 

the interaction term is 

V (x) = (1 - x) 2 v n (0) + 2x(l - x)v 12 (0) + x 2 v 22 (0), (31) 

and the contribution due to the matrix is 

Q{x) = (1 - xHo(O) + awao(O). (32) 

The latter has a particularly simple form; it is independent of p b and only linearly dependent 
on x. This leads to several special features of the fluid inside the matrix. First, the pressure 
of the annealed GCM fluid in the matrix is independent of the matrix density p b . As the 
pressure P is obtained from 

where v = l/p b is the volume per particle and the matrix contribution to the free energy 
(jinj) is independent of p b , it follows that P is independent of p b Q . 

Second, the phase boundaries in the (x, p b ) plane of the fluid inside the matrix are the 
same as those without the matrix (it is known that the binary GCM fluid can exhibit fluid- 

□ n 

fluid phase separation 20, 21]). This can be seen as follows: the locus of the spinodal line 
can be obtained from the Gibbs free energy per particle, gu = fn + Pv, via the condition: 

'd 2 g 12 



dx 



2 



= 0. (34) 

P,T 



Since P is independent of p h and the terms in /12 (and therefore in #12) that are functions of 
Pq are only linearly dependent on x, we see from (|34j) that in the (x, p b ) plane the spinodal 
line will be the same for all values of pg. The binodal (coexistence curve) is obtained from 
the conditions of equal pressure P and equal chemical potentials in the coexisting phases, 

P a = Pp 

P^l,a — pl,/3 

P2,a = P2,/3, (35) 

where a and (3 label the coexisting phases. From Eq. ()27|) . we see that the presence of the 
matrix results in only a linear shift in the chemical potentials of the fluid particles, and this 
is the same in both fluid phases, a and (3. Therefore, the binodal curve for the binary fluid 
in the matrix is at the same location in the (x, p b ) plane as it is without the matrix. 

The special form of the free energy makes for a simple analysis of the situation when 
one has a bulk binary GCM fluid acting as a reservoir and coupled to the system inside the 
matrix that consists of quenched particles of a third species of GCM particles. The chemical 
potentials, pi and P2, are set by the fluid in the bulk reservoir. One can calculate p\ and 
P2, the densities of the fluid in the matrix, using Eq. (j2Zj), i-e. the densities of the fluid in 
the matrix are simply those of a bulk fluid with chemical potentials equal to p^, n = 1,2, 
given by Eq. (12*71) . 

For a binary GCM fluid exhibiting phase separation, this means that the chemical po- 
tentials, correspond to a fluid on one side of the binodal of the reservoir (say, 
the a-phase), but the effective chemical potentials , given by Eq. (J2"7)) . for the fluid 
inside the matrix can correspond to the other side of the binodal (the /3-phase). Under 
these circumstances a (capillary) condensation of the coexisting /3-phase occurs in the ma- 
trix. Such condensation transitions are well-known for fluids in porous media (the matrix) 
DDDDDD Perhaps the best known case is that of a gas condensing to a liquid at 
pressures below the bulk saturated vapour pressure or, for fixed pressure, at temperatures 
above the bulk boiling temperature. In this case the bulk a-phase corresponds to the gas 
and the /3-phase to the liquid j^ . 

Owing to the simplicity of the present model fluid and the theoretical description that 
we use, Eq. (J2*T|) provides a very simple mapping between points in the (x, p b ) plane corre- 
sponding to the density and concentration of the fluid in the bulk reservoir and the same 
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FIG. 1: The bulk phase diagram for a binary mixture of GCM particles with eyij^w = 0.944 and 



R22/R11 = 0.665 (see also Ref. 21]) in terms of the total density, p b , and the concentration, x, 
of the smaller species 2. In a) the paths A r and B r (arrows) are at constant total density in the 
reservoir. The paths A m and B m indicate the corresponding paths for the fluid in a matrix of 
density p^Rfi = 10~ 3 . b) shows an expanded view of the phase diagram. The three additional 
pairs of binodal (solid line) and spinodal (dashed) lines show the loci in the reservoir phase diagram 
which correspond [via Eq. (|27|) ] to the binodal/spinodal of the fluid adsorbed in the matrix, plotted 
for three different matrix densities. In a) and b) the solid line whose ends are denoted by filled 
circles is the thin-thick adsorbed film transition of the binary fluid adsorbed around a single big 
GCM particle with pair potenti al p arameters eoi = &bT, £02 = 0.8/cbT, Rqi/Ru = 5.0 and 
R02/R11 = 4.97 - see Refs. [23] • The two points (o) are where the spinodal for this single 
particle thin-thick adsorbed transition meets the bulk fluid binodal. 



quantities for the fluid in the matrix. For a given matrix density (and given parameters of 
the pair potential between matrix and fluid particles) we can easily determine the matrix 
condensation phase behaviour. 



V. RESULTS FOR A PARTICULAR BINARY MIXTURE 

We consider cases where the annealed fluid components are characterised by the pair 
potential parameters en = e 2 2 = 2ksT, €12 = 1.8877/cbT, R22 = 0.665-Rn; this implies from 
Eq. (jH)) that Ru = 0.849-Rii- This particular binary mixture exhibits bulk fluid-fluid phase 
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separation and was investigated in a number of previous studies |2jJ,|2i 12^12^0- The 
bulk phase diagram is displayed in Fig. ^ The parameters chosen for the pair potentials 
characterising matrix particles are e o = 2k B T, e i = k B T, e 02 = 0.8k B T, R 00 = 7Rn, 
from Eq. (|21j) it follows that i? i — 5i?n and Rq2 = 4.97-Rn. The reason for this choice of 
parameters is that much is known about the behaviour of an annealed fluid of such (big) 



particles dissolved in the binary (solvent) fluid |25|, |26j, |2jj ; hence we can make connections 
with that body of work. In Fig.^b) we display a portion of the phase diagram, together with 
the loci in the reservoir phase diagram of the state points which correspond [via Eq. ([27])] to 
points on the binodal/spinodal for the fluid adsorbed in the matrix. We display results for 
three densities of the quenched particles, p b Q R\i = 0.0001, 0.001 and 0.002 (or, equivalently, 
Pq-Roo = 0.0343, 0.343 and 0.686). As the matrix density p b is increased, the binodal for the 
fluid adsorbed in the matrix shifts to larger x, away from the reservoir binodal. Note that 
one would obtain the same shift for a matrix with a different set of pair potential parameters 
and matrix densities, so long as the products Po^on(O) remain the same (recall Eq. flZZj) )- 
For the above set of matrix pair potential parameters we know 

y, 

a single big 

particle can be surrounded by a thick adsorbed film of the coexisting phase. Furthermore, 
on decreasing x, for certain fixed total fluid densities, there can be a transition from a thin to 
a thick adsorbed film, see Fig. ^ Given this information, and the typical size of the (thick) 
adsorbed film, one would expect the matrix condensation line (e.g. those in Fig. Ub)) to be 
located in the vicinity of the single particle thin-thick adsorbed film transition line, provided 
the density of the quenched matrix particles is sufficiently high that their adsorbed films 



(start to) overlap. I s 
a radius ~ 7Rn t 



ear bulk coexistence, the adsorbed film around a single big particle has 



26 



271 ]. One would therefore expect the matrix condensation line to 



be near the single particle thin-thick adsorbed film transition line when the matrix density 
Pq ~ (47r(7i?n) 3 /3) _1 = 0.7 x lO -3 /?^ 3 . This is indeed what we observe. In fact, we find 
that the spinodal line for the single particle thin-thick transition line lies very close to the 
spinodal line (in the reservoir phase diagram) for a matrix of density p b Q = lO -3 ^?^ 3 - see Fig. 
Eb- This observation provides physical insight into the mechanism driving the condensation 
in the matrix: condensation occurs when the big matrix particles are sufficiently close that 
the adsorbed films around each of them can overlap. This means that condensation in the 
present soft core system is somewhat different from the condensation that we alluded to 
earlier where the condensation occurs in the gaps in the matrix 1.6. cLS cL surface effect, 
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FIG. 2: Fluid densities, p w inside the matrix corresponding to a reservoir with fixed total density 
p b R\i = 8.0 and concentration x = 0.96 plotted versus (uniform) matrix density p$. For /Cg-Rfi > 
2.3 x 10~ 3 , condensation of the coexisting fluid phase occurs in the matrix. Note that p\ is a 
non-monotonic function of Pq, whereas p\ and the total density p b , are monotonically decreasing 
functions of p^. 

where the films on (opposite) matrix surfaces join. Rather, for the present soft-core system, 
the onset of condensation is related to the density profile of of the solvent around and inside 
a single matrix particle. One can therefore determine much about the fluid behaviour in the 
bulk matrix from knowing the behaviour around a single matrix particle. 

In Fig. E]we display the fluid densities, p h n1 n = 1, 2, inside a matrix of quenched particles 
with uniform density for the case of a bulk fluid reservoir with total density p b R\i = 8.0 
and concentration x = 0.96. Results are plotted for a range of matrix densities p h Q . For 
p\)R\i > 2.3 x 10~ 3 , we find condensation of the coexisting fluid phase in the matrix. We 
also find that p\ (dashed-dotted line) is a non-monotonic function of p$, whereas p\ (dashed 
line) and the total density p b = p\ + p\ (solid line), are monotonically decreasing functions 
of Po- F° r ^0.1 the fluid densities are very small; both species are expelled from the 

matrix. This behaviour is typical of that found for reservoir state points lying on the right 
hand side of the bulk binodal in Fig. ^ 
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VI. INHOMOGENEOUS FLUID DENSITY PROFILES: PLANAR INTERFACE 



In this section we determine the inhomogeneous fluid density profiles of the binary fluid 
for a simple model of an interface in which the quenched matrix particles are described by 
an average density distribution: 



i.e. the matrix particles are confined to the half-space with (Cartesian) coordinate z < 
and the fluid reservoir is at z > 0. The planar interface is located at z = 0, where z is the 
coordinate perpendicular to the interface. The step function in Eq. (JHUj) is not the most 
realistic model for the matrix density profile; one would expect some structure in the profile 
for z — ► that will depend on how the quenched matrix is formed. However, the simple 
choice (XM)\ allows us to explore some of the interfacial phenomena that can be exhibited by 
the present system. We do not expect our results to be changed qualitatively were we to 
chose a more realistic profile. 

In Fig. El we display the density profiles p n (z), n = 1, 2, of the binary GCM fluid for cases 
where the reservoir fluid has a fixed bulk density p b R\i = 7.0 and concentration x = 0.888 
that is near to coexistence (located at x coex = 0.88544). The density profiles are calculated 
for six matrix densities Pq-^ii = 10~ 5 , 10~ 4 , 10~ 3 , 10~ 2 , 0.05 and 0.1. Provided the matrix 
density is sufficiently high the coexisting fluid condenses in the matrix. For large values of 
the matrix density the fluid is expelled from the bulk of the matrix - observe the results for 
PqR^i = 0.1. The densities of the fluid in the bulk of the matrix (z — ► — oo) behave in a 
similar manner to the results displayed in Fig. El Strikingly, in addition to condensation in 
the bulk of the matrix, a wetting film of the coexisting phase develops on the surface of the 
matrix. This film is present in spite of the fact that the fluid density in the matrix is very 
low, owing to the high matrix density. Similar wetting behaviour was found for the present 
GCM fluid mixture at a planar wall, for several choices of (repulsive) wall-fluid potentials 



Since the decay of the effective matrix potential ()28j) into the bulk of the reservoir fluid 
is short ranged, the thickness I of the wetting film of the coexisting phase on the surface of 
the matrix (i.e. I is the distance of the wetting film interface from z = 0), should increase as 




(36) 
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FIG. 3: Density profiles of the binary GCM fluid at the planar interface between a matrix with 
density profile given by Eq. (|36|) and a reservoir of fluid (at z > 0) with fixed total bulk density 
p b Rfi = 7.0 and concentration x = 0.888. The density profiles for species 1 are shown in a), and 
those for species 2 in b). The density profiles are calculated for matrix bulk densities Po-^-ii = 10 -5 , 
10 -4 , 10 -3 , 10 -2 , 0.05 and 0.1. Note that for intermediate densities the coexisting fluid phase is 
condensed in the matrix. At this reservoir state point, there is a wetting film extending from the 
interface at z = into the bulk reservoir fluid for the higher matrix densities. 

I ~ — £ w \n(5x), where 5x = x — x coex is the difference in concentration between the reservoir 
concentration x and the concentration at bulk coexistence x coex . £ w is the true correlation 
en^th in the bulk coexisting phase which is wetting the wall (matrix)-reservoir interface 

We investigate the wetting behaviour in more detail in Figs. 0] and 03 where we display 
the density profiles of both species for the binary fluid at the planar interface given by Eq. 
(fHfjj) with fixed p h Q R\i = 10~ 3 . Results are presented for reservoir state points corresponding 
to points along paths A r and B r in Fig. ^a). In Fig. 0] the profiles are calculated along path 
A r (p b Rl 1 = 7) for a range of concentrations approaching the reservoir binodal. Between 
x = 0.937 and 0.934 the coexisting phase condenses in the matrix - note the jump in the 
profiles. As x is decreased further, a wetting film develops on the surface of the matrix. 
This increases in extent as x approaches its value at coexistence. In Fig. El we display the 
density profiles for the same matrix density but now for reservoir total density p b R\i = 6.0, 
corresponding to path B r in Fig. |T]a). Once again the profiles are calculated for a range of 
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FIG. 4: Density profiles of the binary fluid at the planar interface between a matrix with density 
profile given by Eq. (|36j) . with Po-Rfi = 10~ 3 . The bulk fluid reservoir (z — > oo) has fixed total 
density p b R\i = 7.0. The profiles are calculated for concentrations: x = 0.99, 0.96, 0.94, 0.937, 
0.934, 0.93, 0.91, 0.89, 0.886 and 0.88545 (bulk coexistence is at x coex = 0.88544). These correspond 
to a series of points along path A r in Fig.^a). Note that between x = 0.937 and 0.934, condensation 
of the coexisting phase occurs in the matrix. Decreasing x further, a wetting film develops on the 
surface of the matrix. 

concentrations approaching the reservoir binodal and we observe a wetting film developing 
on the surface of the matrix. Note the development of a shoulder in the profile for x = 0.785, 
close to coexistence. Contrary to the case in Fig. 01 the density of the fluid adsorbed in the 
matrix changes continuously as x is varied; there is no jump. In both cases the wetting 
film thickens as x — > x coex , and we find that the film thickness I ~ — £ w ln(5x). In general, 
when the coexisting phase is condensed in the matrix and the reservoir state point is near 
to coexistence we expect the matrix-reservoir interface to be wet completely by a film of 
the coexisting phase asi-> x coex . 

VII. CONCLUSIONS 

We have considered the properties of a ternary GCM mixture, one of whose species is 
quenched, modelling a rigid matrix of (crosslinked) polymers. The two annealed species then 
model a binary polymer solution that is immersed in the matrix. Using a DFT treatment 
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FIG. 5: Density profiles of the binary fluid at the planar interface between a matrix with density 
profile given by Eq. (JHSJ, with p^R^ = 1(T 3 . The bulk fluid reservoir (z — > oo) has total density 
p b R 3 n = 6.0. The profiles are calculated for concentrations x = 0.99, 0.95, 0.92, 0.9, 0.89, 0.87, 0.84, 
0.81, 0.8, 0.79, 0.788, 0.786 and 0.785 (bulk coexistence is at x coex = 0.7848). These correspond 
to a series of points along path B r in Fig. ^ a )- For decreasing x, a wetting film develops on the 
surface of the matrix. Note that, contrary to the case in Fig. the density of the fluid adsorbed 
in the matrix changes continuously with x. 

we have calculated both the behaviour of the binary mixture in the bulk matrix and its 
properties at a planar interface, described by the profile (JHrJ|) . between the matrix and a 
reservoir of the fluid. 

In the first case where the fluid is adsorbed in the bulk matrix with uniform density p^, we 
find that in the total density versus concentration (x, p b ) representation of the adsorbed fluid 
phase diagram, the location of the fluid-fluid binodal is independent of p b Q . This somewhat 
surprising result is a direct consequence of the soft core nature of both the adsorbed fluid 
and matrix particles as described within a second-order (in densities) free energy - see Eq. 
f!29|) . In particular the contribution due to the matrix is independent of the density of the 
annealed fluid. For given fluid densities in the reservoir, the calculation of the fluid densities 
in the matrix is very straightforward, and mapping out lines of (capillary) condensation in 
the matrix, in the (x, p b ) phase diagram, is readily performed. 

Applying a Helmholtz free energy functional more accurate than the simple RPA func- 
tional we would expect the fluid phase diagram in the (x, p b ) plane to depend explicitly 
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on the matrix densit y £>n • However, bearing in mind the accuracy of (|23j) at high fluid den- 
sities 19, 0, H Q Q, 

we expect only a weak dependence on p b Q in the true phase 



diagram of the binary GCM immersed in a GCM matrix. 

In earlier work [2|J |2y, |22| , we found that due to the soft-core nature of the particles, 
condensation of the coexisting phase occurs both on the surface and inside a single big 
particle. Making connection with this work, we found that for matrices composed of big 
GCM particles with density Po-Roo ~ 1 0- e - densities near the overlap concentration), the 
mechanism for capillary condensation of the coexisting fluid phase in the matrix is the joining 
of the thick films of the coexistin g ph ase that form around individual big GCM particles 



as bulk coexistence is approached [2jj |2y, |27j . This mechanism is somewhat different from 
what occurs in hard core systems, where the wetting films grow on the surface of big matrix 
particles and condensation occurs throughout the matrix by the joining of these surface 
films. 

Much of the simplicity inherent in the present theory stems from the fact that the effect of 
the matrix on the fluid can be described in terms of a shift in the effective chemical potential 
[c.f. Eq. (12 7j) ] and an effective one body external potential - see Eq. (|28jl . Note that this 
simplification is not clear at the outset: a particular configuration of matrix particles can, of 
course, be treated as an external potential for the adsorbed fluid, but having averaged over 
the ensemble of possible matrix configurations, the description of the matrix then enters on 
the level of the (ensemble average) one body density profile of the matrix particles. That 
this problem can be subsequently mapped to an effective one body potential, which is much 
simpler to treat than the potential exerted by a particular configuration of matrix particles, 
makes determining properties of the inhomogeneous fluid in contact with an inhomogeneous 
matrix tractable. Following this procedure we were able to calculate adsorption and wetting 
behaviour of the binary mixture on the surface of the matrix in contact with a reservoir of the 
fluid. In the case of a reservoir at state points near to coexistence, we found a rich interplay 
between the fluid condensation in the bulk of the matrix and wetting behaviour at the 
interface. We conclude that the present DFT approach constitutes a powerful approach for 
treating fluids in inhomogeneous random media. Although the GCM model is a very simple 
one, it is physically realistic and the results obtained from analysing its properties should 
provide much insight into bulk and interfacial phenomena of fluids adsorbed in random 
porous media. 
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